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Abstract. Recent analyses of several isospin effects in heavy- ion reactions have 
allowed us to constrain the density dependence of nuclear symmetry energy at 
ly-^ I sub-saturation densities within a narrow range. Combined with constraints on 

£f) . the Equation of State (EOS) of symmetric nuclear matter obtained previously 

' from analyzing the elliptic flow in relativistic heavy-ion collisions, the EOS of 

\Q . neutron-rich nuclear matter is thus partially constrained. Here we report effects 

' of the partially constrained EOS of neutron-rich nuclear matter on the mass- 

radius correlation, moment of inertia, elliptical deformation and gravitational 
radiation of (rapidly) rotating neutron stars. 
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1. EOS of neutron-rich nuclear matter partially constrained 
by heavy-ion reactions 

The EOS of isospin asymmetric nuclear matter can be written within the well-known 
parabolic approximation as 



E(p, S) = E(p, 5 = 0) + E sym (p)S 2 + 0(S 4 ), 



(1) 
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where 8 = (p n — p P )/(pp + p n ) is the isospin asymmetry with p n and p p denoting, 
respectively, the neutron and proton densities, E(p, S — 0) is the EOS of symmetric 
nuclear matter, and E sym (p) is the density-dependent nuclear symmetry energy. 
The latter is very important for many interesting astrophysical problems [[H [2], the 
structure of rare isotopes [[3] and heavy-ion reactions [ [H [5] [6] [7l [8] . However, the 
density dependence of the nuclear symmetry energy has been the most uncertain 
part of the EOS of neutron-rich matter. Fortunately, comprehensive analyses of 
several isospin effects including the isospin diffusion [TUj and isoscaling [[TT] in 
heavy- ion reactions and the size of neutron skin in heavy nuclei [|T5] have allowed us 
to constrain the density dependence of the symmetry energy at sub-saturation densi- 
ties within approximately 31.6(p/po) ' 69 an d 31.6(p/po) 105 a s labelled by x = and 
x = —1, respectively, in the lower panel of FigfT] [ [T3l [14] . While these constraints 
are only valid for sub-saturation densities and still suffer from some uncertainties, 
compared to the early situation they represent a significant progress in the field. 
Further progress is expected from both the parity violating electron scattering ex- 
periments [ 15 at the Jefferson lab that will help pin down the low density part of 
the symmetry energy and heavy-ion reactions with high energy radioactive beams at 
several facilities that will help constrain the high density behavior of the symmetry 
energy [|S]. 

For many astrophysical studies, the EOS is usually expressed in terms of the 
pressure as a function of density and isospin asymmetry. Shown in Fig. [T] are 
the pressures for two extreme cases: symmetric (upper panel) and pure neutron 
matter (lower panel). The green area in the density range of 2 — 4.6po is the ex- 
perimental constraint on the pressure Pq of symmetric nuclear matter extracted 
by Danielewicz, Lacey and Lynch from analyzing the collective flow data from 
relativistic heavy- ion collisions [ [6] . It is seen that results from mean- field calcula- 
tions using the phenomenological momentum-dependent (MDI) interaction [[16] , the 
Dirac-Brueckncr-Hartree-Fock approach with the Bonn B potential (DBHF) [I17j. 
and the variational calculations by Akmal, Pandharipande, and Ravenhall (APR)[ 
118] are all consistent with this constraint. For pure neutron matter, its pressure 
is Ppnm = Pq + P 2 dE sym / dp and depends on the density dependence of nuclear 
symmetry energy. Since the constraints on the symmetry energy from terrestrial 
laboratory experiments are only available for densities less than about 1.2po as in- 
dicated by the green and red squares in the lower panel, which is in contrast to the 
constraint on the EOS of symmetry nuclear matter that is only available at much 
higher densities, the most reliable estimate of the EOS of neutron-rich matter can 
thus be obtained by extrapolating the underlying model EOS for symmetric matter 
and the symmetry energy in their respective density ranges to all densities. Shown 
by the shaded black area in the lower panel is the resulting best estimate of the 
pressure of high density pure neutron matter based on the predictions from the MDI 
interaction with x=0 and x=-l as the lower and upper bounds on the symmetry 
energy and the How-constrained EOS of symmetric nuclear matter. As one expects 
and consistent with the estimate in Ref . [ [6] , the estimated error bars of the high 
density pure neutron matter EOS is much wider than the uncertainty range of the 
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Fig. 1. (Color online) Pressure as a function of density for symmetric (upper panel) 
and pure neutron (lower panel) matter. The green area in the upper panel is the 
experimental constraint on symmetric matter. The corresponding constraint on 
the pressure of pure neutron matter obtained by combining the flow data and an 
extrapolation of the symmetry energy functionals constrained below 1.2po by the 
isospin diffusion data is the shaded black area in the lower panel. Results taken 
from Refs. [© Q3]. 
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EOS of symmetric nuclear matter. For the four interactions indicated in the figure, 
their predicted EOS's cannot be distinguished by the estimated constraint on the 
high density pure neutron matter. In the following, the astrophysical consequences 
of this partially constrained EOS of neutron-rich matter on the mass-radius corre- 
lation, moment of inertia, and the elliptical deformation and gravitational radiation 
of (rapidly) rotating neutron stars are briefly discussed. More details of our studies 
on these topics can be found in Refs. [fT9l I20 t I2T ] 122 ] [23] . 



2. Nuclear constraints on the mass-radius correlation, moment 
of inertia, elliptical deformation and gravitational radiation 
of rapidly rotating neutron stars 
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Fig. 2. (Color online) Gravitational mass versus equatorial radius for neutron stars 
rotating at v = 716 Hz and f = 1122 Hz. Taken from Rcf. [|2"2]. 

The partially constrained EOS of neutron-rich nuclear matter has important 
ramifications on properties of neutron stars. As a first example, in Fig. [5] we show 
the mass-radius correlations for the two fastest rotating neutron stars known as of 
today. These pulsars spin at 716 [[24] and 1122 Hz [[25], respectively. However, 
based only on the observational data available so far, their properties have not 
yet been fully understood. The analysis of their properties based on the EOS and 
symmetry energy constrained by the terrestrial laboratory data is thus especially 
interesting. Setting the observed frequency of the pulsar as the Kepler frequency, 
corresponding to the highest possible frequency for a star before its starts to shed 
mass at the equator, one can obtain an estimate of its maximum radius as a function 
of mass M, 

( M \ 1/3 

iWAf ) = X t^tt- km, (2) 
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Fig. 3. Gravitational-wave strain amplitude as a function of the neutron star mass. 
The error bars between the x = and x = — 1 EOSs provide a limit on the strain 
amplitude of the gravitational waves to be expected from these neutron stars, and 
show a specific case for stellar models of IAMq. Taken from ref . [ fT9] . 



with x — 20.94 for rotational frequency v — 716 Hz and x = 15.52 for v = 1122 Hz. 
The maximum radii are shown with the dotted lines in Fig. [2] It is seen that the 
range of allowed masses supported by a given EOS for rapidly rotating neutron 
stars becomes narrower than the one of static configurations. This effect becomes 
stronger with increasing frequency and depends upon the EOS. Since predictions 
from the x = and x = —I EOSs represent the limits of the neutron star models 
consistent with the experimental data from terrestrial nuclear laboratories, one can 
predict that the mass of the neutron star rotating at 1122 Hz is between 1.7 and 
2.1 solar mass [I22j. 

Another interesting example is the gravitational radiation expected from el- 
liptically deformed pulsars. Gravitational waves (GWs) are tiny disturbances in 
space-time and are a fundamental, although not yet directly confirmed, prediction 
of General Relativity. Gravitational wave astrophysics would open an entirely new 
non-electromagnetic window to the Cosmos, making it possible to probe physics 
that is hidden or dark to current electromagnetic observations [ 126] . Elliptically de- 
formed pulsars are among the primary possible sources of the GWs. Very recently 
the LIGO and GEO collaborations have set upper limits on the GWs expected from 
78 radio pulasrs [ [27] . Gravitational waves are characterized by a strain amplitude 
ho which can be written as 



with x = -\/20487r 5 / 15G / c 4 . In the above equation, r is the distance between 
the pulsar and the detector, and the $22 is the quadrupole moment of the mass 
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distribution. For slowly rotating neutron stars, one has [ 128) 

— ^(^(^H^)" (4) 

In the above expression, a is the breaking strain of the neutron star crust which 
is rather uncertain at present time and lies in the wide range a = [10 -5 — 10 -2 ] [ 
129] . In our estimate, we use the maximum breaking strength, i.e. a = 1CP 2 . In 
Fig. [3]we display the GW strain amplitude, ho, as a function of stellar mass for three 
selected millisecond pulsars which are relatively close to Earth (r < OAkpc) and have 
rotational frequencies below 300 Hz. It is interesting to note that the predicted ho 
is above the design sensitivity of LIGO detector. The error bars in Fig. [3] between 
the x = and x = — 1 EOSs provide a constraint on the maximal strain amplitude 
of the gravitational waves emitted by the millisecond pulsars considered here. The 
specific case shown in the figure is for neutron star models of \AMq. Depending 
on the exact rotational frequency, distance to detector, and details of the EOS, the 
maximal ho is in the range ~ [0.4 — 1.5] x 10~ 24 . These estimates do not take 
into account the uncertainties in the distance measurements. They also should be 
regarded as upper limits since the quadrupolc moment (Eq. ([?])) has been calculated 
with a = 10~ 2 (where a can go as low as 10~ 5 ). 

To emit GWs a pulsar must have a quadrupole deformation. The latter is nor- 
mally characterized by the ellipticity which is related to the neutron star maximum 
quadrupole moment $22 and the moment of inertia via [ [28] 




(5) 

For slowly rotating neutron stars, one can use the following empirical relation [ 130] 



(0.237 ±0.008)Mi? 2 



, n Mkm ( Mkm^ ' 
1 + 4.2 h90 

M G R \M q R 



(6) 



This expression is shown to hold for a wide class of equations of state which do 
not exhibit considerable softening and for neutron star models with masses above 
1M Q [[30]. 

Fig. 0] displays the neutron star moment of inertia (left panel) and ellipticity 
(right panel). It is interesting to mention that a fuducial value of I zz — 10 45 g cm 2 
is normally assumed in the literature. Our calculations indicate that the I zz is 
strongly mass dependent. This observation is consistent with previous calculations. 
Moreover, the ellipticity decreases with increasing mass. The magnitude is above the 
lowest upper limit of 4 x 10~ 7 estimated for the PSR J2124-3358 [ 27J. Interestingly, 
essentially all obervables depend strongly on the EOS of neutron-rich matter. In 
particular, the MDI EOSs, adopting the same symmetric matter EOS but different 
density dependence of the symmetry energy, sets useful nuclear boundaries for these 
gravitational wave observables. 
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In summary, the heavy-ion physics community has made significant progress in 
constraining the EOS of neutron-rich nuclear matter in recent years. In particular, 
comprehensive analyses of several isospin effects including the isospin diffusion and 
isoscaling in heavy-ion reactions and the size of neutron skins in heavy nuclei have 
allowed us to constrain the symmetry energy at sub-saturation densities within ap- 
proximately 31.6(p/po) ' 69 and 31.6(p/po) 105 - While the currently existing data 
only allowed us to constrain the symmetry energy and thus the EOS of neutron- 
rich matter in a narrow range, it can already help to put some useful constraints 
on several interesting observables in astrophysics, such as the mass-radius correla- 
tion, moment of inertia, and the elliptical deformation and gravitational radiation 
of (rapidly) rotating neutron stars. With the parity violating electron scattering 
experiments and heavy-ion reactions with high energy radioactive beams, it will be 
possible in the future to map out accurately the entire density dependence of the 
symmetry energy. 
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